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We summarize our current understanding of instantons in the large N c limit of 
QCD. We also present some recent results from simulations of the instanton liquid 
in QCD for N c > 3. 



1 Historical Overview 

Both instantons and the large N c expansion are important corner stones in 
our attempts to understand non-perturbative aspects of QCD. Instantons are 
related to topological aspects of QCD and play and important role in under- 
standing the U(\)a anomaly. There is also a growing body of evidence that 
the mechanism for chiral symmetry breaking in QCD is intimately connected 
with instantons. The large N c expansion is based on the idea that 1/N C can 
be used as expansion parameter in QCD. The main assumption is that QCD 
with many colors is a theory that exhibits confinement and chiral symmetry 
breaking, and that this theory is smoothly connected to the real world, in 
which N c is equal to three. There is an impressive amount of evidence, mainly 
of phenomenological nature, but also from the lattice, that this assumption is 
indeed correct. 

It is commonly believed, however, that instantons and the large N c expan- 
sion are fundamentally incompatible with each other. This view is maybe best 
explained by providing a short review of the literature on the subject. The in- 
stanton solution was discovered^ in 1975. In the following year 't Hooft found 
the fermion zero mode in the background field of an instanton and explained 
the connection between instantons and the U(1)a anomalyEI. In 1978 Witten 
wrote a very influential paper U entitled "Instantons, the Quark Model, and 
the 1 /N expansion" . He noted that straightforward N c counting suggests that 
the rj' mass squared scales as 1/N C , whereas classical effects, such as instan- 
tons, scale as exp(— 1/g 2 ) ~ exp(— N c ). He explicitly stated that in trying to 
understand the r/' mass and other non-perturbative phenomena in QCD "it is 
necessaxy to choose between the large N expansion and instantons" . In 1979 
WittenQ and Venezianou derived a relation between the mass of the r( and the 
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topological susceptibility in pure gauge theory, 



2N, 



(m„ + m n , - 2m K ) = xtop- (1) 



Witten suggested that for all N c instantons "melt" and cannot be distinguished 
from perturbative fluctuations. He also proposed that the topological suscep- 
tibility is of order O(l) in the large N c limit, and that it is dominated by 
non-perturbative fluctuations other than instantons, presumably related to 
confinement. Indeed, Veneziano's paper is entitled "£7(1) without Instantons". 

Simple N c counting implies that both the t/' — 7r splitting and the p — uj 
splitting are of order 1/N C . In the real world, of course, these two phenomena 
are very different in magnitude. This suggests that there is some physical 
effect in the rf channel that leads to unusually large 1/N C corrections. This 
observation was discussed in more detail by NovikovB et al. in a paper entitled 
"Are all hadrons alike?" . Novikov et al. observed that large corrections to the 
OZ1 rule appear many channels, not just pseudoscalar mesons, but also scalar 
mesons and glueballs. They also noticed that these are precisely the same 
channels in which direct instanton effects appear. 

This observation led to the formulation of the instanton liquid model of the 
QCD vacuumQtrB. The instanton liquid model postulates that the N c = 3 QCD 
vacuum is populated by localized, approximately dual or anti-self dual lumps, 
instantons. The density of instantons is approximately (N/V) ~ 1 fm -4 while 
the size is p ~ 1/3 fm. These numbers reproduce the topological susceptibility 
in the pure gauge theory xtop — (200 MeV) 4 and the chiral condensate (V>"0) — 
— (230 MeV) 3 . More detailed calculations show that the instanton liquid model 
successfully .describes an impressive amount of date on hadronic correlation 
functionsBBB. 

In the 1980's researchers also began to study the topological structure of 
QCD on the lattice. It was found that the topological susceptibility in pure 
gauge QCD isli3 xtop — (200 MeV) 4 , as predicted by the Witten- Veneziano rela- 
tion equ. (0). However, it was also observed that the topological susceptibility 
is very stable under cooling, and appears to be dominated by semi-classical 
configurations. Lattice simulations alsn appear to confirm the values of the 
key parameters of the instanton liquidcl, (N/V) ~ 1 fm~ 4 and p ~ 1/3 fm. 

More recently lattice simulations have started to focus on the structure of 
low-lying eigenvectors of the Dirac operator and the mechanism of chiral sym- 
metry breaking. The instanton model predicts that the lowest eigenstates of 
the Dirac operator, which dominate chiral symmetry breaking, are linear com- 
binations of localized, approximately chiral states associated with the fermionic 
zero modes of individual instantons and anti-instantons til. This picture has 
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been confirmed by lattice calculationslHlO, although there is some controversy 
concerning the question whcth£r_thc size of the chiral lumps is in agreement 
with the instanton predictionEj'LD. 

2 Instantons and the OZI rule 

In this section we would like to remind the reader how instantons lead to large 
violations of the OZI rule, 't Hoofn explained that the effect of an instanton 
on fermionic correlation functions can be summarized in terms of an effective 
interaction 

L = G [(^T a ij) 2 - (V^) 2 - (^l5Ta^) 2 + {^1^?} , (2) 

where G is related to the tunneling amplitude and t is an isospin matrix. 
On the single instanton level we can directly read off the interaction in the 
channel characterized by the current ipTip. The interaction is attractive in the 
pion (ipi"/§T?p) and sigma {ipip) channels, and repulsive in the eta prime (ip'il5' l l J ) 
and ao (-ipTip) channels. We note that the OZI violating interaction the rj' and 
a channel is not suppressed with respect to the OZI allowed contribution in 
the 7r and ao sector. 

It is worth repeating why this is so. The instanton interaction corresponds 
to the contribution of fermionic zero modes to the quark propagator. Since 
there is exactly one zero mode for every flavor, there are no diagonal (uu)(uu) 
or (dd)(dd) interactions. Second, since the fermion zero modes for quarks and 
anti-quarks have opposite chirality, the interaction is also off-diagonal in the 
basis spanned by right and left-handed fermions q^, q^- As a consequence 
the sign of the interaction flips in going from the scalar to the pseudoscalar 
channel, and in going from / = to I = 1 states. Also, to leading order there 
is no instanton contribution to the interaction in vector channels. 

In order to be more quantitative we consider correlation functions involving 
the scalar and pseudoscalar currents introduced above. The correlators are 
defined by 

U„(x,y) = (Tr (S(x,y)ij 5 S{y,x)ij 5 )), (3) 
Il s (x,y) = (Tr (S(x,y)S(y,x))}, (4) 
n n '(a;, y) = (Tr (S(x, y)i-f 5 S(y, x)ij 5 )) - 2(Tr (i-/ 5 S(x, x)) Tr (i-f 5 S(y, y)))(5) 
n CT (x, y) = (Tr (S(x, y)S(y, x))) - 2(Tr (S(x, x)) Tr (S(y, y))), (6) 

where S(x,y) is the fermion propagator and (.) denotes the average over all 
gauge configurations. The OZI-violating difference between the ir, r/' and <r, ao 
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Figure 1: Correlation functions in the pion, delta (ao), and rho meson channels. The corre- 
lators are normalized to free field behavior (rio(rr) ~ a; -6 ). 
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Figure 2: OZI violating disconnected correlation functions in the pseudoscalar [rf — 7r), scalar 
(ct — ao), and vector (u> — p) channel. The correlators are normalized as in Fig.l. 
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channels is determined by "disconnected" (or double hairpin) contributions to 
the correlation functions. 

In Figs. 1 and 2 we show results for the correlation functions obtained 
in unquenched instanton simulations of the instanton liquid EJ. We observe 
that the correlation function in the OZI violating rj' — n and a — clq channels 
is as large as the correlation function in the OZI allowed ir channel whereas 
the OZI violating u — p channel is much smaller—This pattern was also found 
in a lattice calculation by Isgur and Thacker It would be interesting to 

perform similar calculations in other channels in which large violations of the 
OZI rule might be present. Examples are the mixing of scalar glueballs with 
scalar mesons or multi-pion states, the AI =1/2 rule, and OZI suppressed 
decays of heavy quark bound states. 



3 Instantons at Large N c 

In the previous section we argued that the OZI violating direct instanton con- 
tribution in the rj' channel is not suppressed compared to the instanton con- 
tribution in the 7r channel. But what is the N c dependence of these effects? 
Witten argued that instanton effects are proportional to exp(— %tt 2 /g 2 ) and 
therefore scale as exp(— N c ). But this is only true for instantons at the cutoff 
scale. In perturbation theory, the tunneling rate is of the form 

dn(p)~p- 5 ( P A) b dp, (7) 

where b — HN c /3 — 2iV//3 is the first coefficient of the beta function. This 
implies that small instantons p < A^ 1 are suppressed as N c — > oo, but large 
instantons are not. 

This is observation is not very useful, because instantons of size A -1 have 
action S ~ 1 and are clearly not semi-classical objects. We can be some- 
what more quantitative, however. The one-loop result contains additional N c 
dependent factors that are related toJjbe collective coordinate measure. Ex- 
ponentiating everything we can write tiJH'Lj 

— ~ exp[N c F(p)], (8) 

with F(p) = 2- S (p) + 21og( S (p)) + . . . and s(p) = S(p)/N c = (87T 2 )/(N c g 2 (p)). 
The function F(p) has a non-trivial zero for s* ~ 5 so the density of instantons 
of size p — p* ~ 0.2A -1 fm is fixed as N c — ■> oo. This is a much more optimistic 
conclusions, because it implies that instantons with action S = N c s* ~ 5N C 3> 
1 can survive in the large N c limit. The precise value of p* and s* depends 
on the renormalization scheme and higher order corrections. A fixed point in 
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the instanton size distribution was indeed observed in the lattice calculation 
of Lucini and TeperEa, see also Fig. 3. However, these authors also find that 
the density of instantons of size p > p* keeps growing as N c —* oo, and that as 
a result there is no well defined average instanton size. 

4 Supersymmetric Theories 

While the fate of instantons in the large N c limit of QCD is still mysterious 
some progress has been made in understanding instantons in the large N c 
limit of supersymmetric gauge theories. The most impressive accomplishment 
is Maldacena's discovery of an explicit master field in the case of N = 4 
supersymmetric QCD at large N c in the limit of large 't Hooft coupling A = 
g N c . The master field is given by supergravity in an AdS$ x S5 background 
together with a set of rules that relate gauge theory to supergravity observables. 
It is very interesting to study how instantons appear in this correspondences. 
On the supergravity side gauge theory instantons appear as D-instantons and 
are characterized by a point on AdS§ x S5. This integration over this point 
corresponds to the collective coordinate integration on the gauge theory side. 
The AdSs part is related to the integration over d 4 xdp/p 5 while the £5 arises 
from the fermion zero modes E3. 

The case of multi-instanton configurations is even more interesting. Naively 
one would think that a fc-instanton contribution involves an integration over 
(AdS$ x S^. In the large N c limit one finds, however, that the saddle 
point configuration is given by k instantons in commuting SU C2j subgroups of 
SU(N C ) all with the same size and located at the same point c3. As a result, 
the collective coordinate measure contains only one copy of AdS§ x S5. The 
multi-instanton configuration is bound by fermion exchanges. 

Progress was also made in understanding N — 2 SUSY QCD. Seiberg and 
Witten determined the low energy effective action of this theory c3. In the 
semi-classical limit, the result can be expressed as the perturbative one-loop 
contribution plus an infinite series on fc-instanton corrections, and, remarkably, 
nothing else. This would seem to imply that the large N c limit of this theory is 
rather boring. Instantons are suppressed as exp(— N c ), monopoles have masses 
0(N C ), and one is left with only the perturbative part. 

This is not correct, however. The methods-.of Seiberg and Witten were 
generalized to arbitrary N c by Klemrn_et al. Ej. This result was analyzed 
in more detail by Douglas and ShenkerEH. The moduli space of the theory is 
characterized by N c — 1 Higgs expectation values. At most points on the moduli 
space the large N c limit is indeed trivial but Douglas and Shenker identified a 
special form of the large N c limit in which instantons and monopoles survive. 
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Figure 3: Instanton size distribution for different values of N c obtained from a numerical 
simulation of the interacting instanton liquid. The instanton size is given in units of the 
Pauli-Vilars scale parameter A ~ 200 MeV. 

5 The Instanton Liquid at Large N c 

In this section we wish to study the question whether it is possible to construct 
a large N c instanton ensemble in QCD that is consistent with standard large 
N c counting. Equ. shows that instantons of action S ~ iV c ^> 1 and size 
p ~ N® can survive in the large N c limit. It is quite natural to assume that the 
total density of these objects would scale as (N/V) ~ N c because instantons 
are essentially SU (2) configurations and the number of mutually commuting 
SU(2) subgroups of SU(N C ) grows as N c . This means that we can have 0(N C ) 
instantons and the instanton liquid remains dilute. Using the trace anomaly 
relation 

(T^) = -^(9 2 G^G- U ) (9) 

with b = UN c /2 — 2Nf/2 we also see that the instanton contribution to the 
vacuum energy scales as £ ~ N£. This is identical to the scaling behavior of 
the perturbative part. 

In order to construct a consistent ensemble of instantons we have to make 
an assumption concerning the fate of large instantons. In the following we shall 
assume that there is a classical 0(Sq ) = 0(N C ) core in the instanton interaction 
for instantons that overlap in group space. This core excludes configurations 
with large or strongly overlapping instantons that are not semi-classical. The 
parameters of the core were fitted to .reproduce the phenomenological values 
of the instanton size and density fort] A^ c = 3. Instanton size distributions 
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Figure 4: Instanton density (N/V), chiral condensate (qq) and topological susceptibility 
obtained from a numerical simulation of the interacting instanton liquid. All quantities are 
given in units of appropriate powers of the Pauli-Vilars scale parameter. 



obtained from numerical simulations of the instanton liquid for iV c = 3, . . . , 6 
are shown in Fig. 3. In these simulations we have assumed that (N/V) ~ 
N c , but this assumption can be verified by computing the free energy of the 
instanton liquid as a function ofE3 (N/V). We note that the distributions 
show the fixed point discussed above and observed in the lattice simulations 
of Lucini and Teper. 

We have also computed the quark condensate and the topological suscep- 
tibility in pure gauge theory, see Fig. 4. We observe that the quark condensate 
scales as (qq) ~ N c . This scaling can be understood using random matrix 
arguments. The quark condensate is approximately given by 

1 /3iV r iV^ 1/2 



w-—,(.-tv) ■ (10) 

The total number of states in the zero-mode zone scales as N zmz ~ (N/V) ~ 

1/2 

N c . This contributes a factor N c to the scaling behavior of (qq). The second 

1 /2 

factor of iV c ' arises from the fact that the quark condensate is inversely pro- 
portional to the average matrix element \Tja\ of the Dirac operator between 
fermion zero modes. Since instanton zero modes live in SU (2) subgroups \Tja | 2 
scales as l/N c . Similar arguments can be used to show that the pion decay 
constant scales as f% ~ N c . 

For a dilute gas of instantons we would expect that the topological sus- 
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ceptibility in pure gauge theory scales as Xtop — (Qto P )/V — ~ N c , 

contrary to the standard large N c assumption \top y-J^c- Further support 
for this assumption was recently provided by Witten Li We have to keep in 
mind, however, that the result Xtop — {N/V) is based on the idea that topolog- 
ical charge fluctuations are Poissonian. Since the classical interaction between 
instantons also grows with N c , this is not necessarily the caseu. Indeed, sim- 
ulations carried out for N c = 3, . . . , 7 seem to indicate that the topological 
susceptibility remains finite in the large N c limit, see Fig. 4. 



6 Instead of Conclusions 

We have argued that instantons provide a very successful explanation of the 
pattern of OZI violation observed in QCD for N c = 3 colors. The fate of 
instantons in the large N c limit remains unclear. It may well be that instantons 
become large and overlap strongly, and the semi-classical description breaks 
down completely. However, it is also possible that the semi-classical picture 
survives. This would seem to require large cancellations and fine tuning, but 
for reasons that we do not understand the parameters of the instanton liquid 
in QCD with N c = 3 colors are already remarkably close to the critical point 
where a smooth large N c limit is possible. 
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